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PULSATION FREQUENCIES AND MODES OF GIANT EXOPLANETS 
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ABSTRACT 

We calculate the eigenfrequencies and eigenfunctions of the acoustic oscillations of giant exoplanets 
and explore the dependence of the characteristic frequency vq and the eigenfrequencies on several 
parameters: the planet mass, the planet radius, the core mass, and the heavy element mass fraction 
in the envelope. We provide the eigenvalues for degree I up to 8 and radial order n up to 12. For the 
selected values of I and n, we find that the pulsation eigenfrequencies depend strongly on the planet 
mass and radius, especially at high frequency. We quantify this dependence through the calculation 
of the characteristic frequency v§ which gives us an estimate of the scale of the eigenvalue spectrum 
at high frequency. For the mass range 0.5 < Mp < 15 Mj, and fixing the planet radius to the Jovian 
value, we find that vq ~ 164.0 x (Mp/Mj) 0AS [iHz, where Mp is the planet mass and Mj is Jupiter's 
mass. For the radius range from 0.9 to 2.0 Rj, and fixing the planet's mass to the Jovian value, 
we find that vq ~ 164.0 x (Rp / 'Rj)~ 2 ' 09 [iR 2, where Rp is the planet radius and Rj is Jupiter's 
radius. We explore the influence of the presence of a dense core on the pulsation frequencies and 
on the characteristic frequency of giant exoplanets. We find that the presence of heavy elements 
in the envelope affects the eigenvalue distribution in ways similar to the presence of a dense core. 
Additionally, we appl y our formalism to J upiter and Saturn and find results consistent with both the 
observational data of lGaulme et al.l ([201 ID and previous theoretical work. 



1. INTRODUCTION 

Pulsation frequencies and modes are potentially useful 
tools with which to study the interior structure of the 
giant planets. Three types of modes are distinguished: 
g-modes are standing internal gravity waves, p-modes are 
standing acoustic waves, and f-modes are of intermedi- 
ate frequency and can be regarded as the fundamental 
mode of either the p- or the g-modes. The correspond- 
ing frequencies are characterized by their radial order 
n and degree I. Several forcing mechanisms may excite 
these modes, maki ng them, in principle , acce s sible to 
direc t observation ([Bercovici fc SchubertI fT987t IMarlevI 
1990). However, the surface movements associated with 
such pulsations are hard to detect and, as of July 2012, 
only Jupite r's global oscillations have been observed . 
The work oflSchmider et all (1199111 IMosser et al.l (119911). 
iMosser et all (|1993l ). and IMosser et all (|2000ft resulted in 
a mea surement of the me an frequency spa cing of 142 ± 3 
uHz (IMosser et al.l 120001 ) . More recently. IGaulme et all 
(|20T1 . using the SYMPA Fourier spectro-imager, de- 
tected Jupiter's global modes with a mean noise level 
five times lower than previously achieved. Their obser- 
vations suggest a mean spacing of 155.3 ± 2.2 fj,Hz. 

Given current technological limitations to the observa- 
tion of giant planet's oscillations, it is reasonable to think 
that the global oscillations of giant exoplanets will not be 
detected in the foreseeable future. However, a theoretical 
exploration of the systematics of the pulsation frequen- 
cies for the broad spectrum of recently discovered giant 
exoplanets might stimulate observers to design methods 
to detect giant planet oscillations, since such oscillations 
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are so diagnostic of structure. 

In this spirit, we calculate the eigenfrequencies and 
eigenfunctions of the pulsationnal modes of giant ex- 
oplanets. We quantify the dependence of the modal 
eigenfrequencies on the planet mass and radius. In addi- 
tion, we focus on the influence of a dense core on these 
quantities, since its presence has al ready been suggested 
in specific extrasola r giant planets (|Burrows et all 120071 ; 
iGuillot et all 12006). Furthermore, we calculate corre- 
sponding models for Jupiter and Saturn themselves, and 
c ompare them wit h pr evious work. 

iVororitsov et all (I1976D . iVorontsov fc Zarkhovl 
(|1981 al ibi) and IVorontsovl (|198ll ) added rotation, 
differential rotation, and ellipticity to their initial 
spherically-symmetric, nonrotating Jovian models. The 
influence of the troposphere on the high-frequenc y 
oscillations was firs t adressed bylVorontsoy et al.l »1 989'1 . 
then in detail by IMosser et al.l (|1994D ~ iProvost et all 
(1993) developed an asymptotic method to determine 
the eigenfrequencies which included the discontinu- 
ity of a Jovian core. They introduced the mean 
spacing or characteristic frequency i/q, defined by: 



2 f ^ 

J planet cq 



where cq is the speed of sound, 

and emphasized the sensitivity of the Jovian oscillation 
spectrum to the presence of a dense core. Since then, 
the Jovian characteristic frequency has been estimated 
to be between 152 and 160 uHz (IProvost et al.l 119931 : 
IGudkovaet aT1ll995l: IGudkova fc Zarkhovlll999ll These 
estimates are consis tent with the recent observations of 
IGaulme et all ([21HTI) . 

Sat urn's os c illatio ns have also been studied theoreti- 
cally. Mar ley (1990) suggested that the f-modes of Sat- 
urn are the most likely to be detected through their po- 
tential influence on that p lanet's rings. U s ing th e tech- 
niques applied to Jupiter, Gud kova et al.l (|1995|) calcu- 
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lated the eigenfrequencies of the lowest-order p-modes of 
Saturn, along with their characteristic frequencies. The 
latter were found to be between 106 and 109 [iHz. 

Throughout our paper, we highlight the general sys- 
tematics and do not take into account the effects o f ro- 
tatio n or oblateness (jVorontsov fc Zarkhovlll98i~aH b1; ILed 
[1991 . Since the a diabatic appro ximation is appropriate 
for Jovian planets (Marlcy 1990), adiabaticity is here as- 
sumed for Jupiter, Saturn, and the entire set of exoplan- 
ets. Finally, since we focus on the giant exoplanet regime, 
we sel ect the appropriate ra nge of measured giant planet 
r adii (lUdrv fc Santosll2007 l) : 0.8 R.j < R P < 2.1 R 7 . 

iGaulme et all (|2011l ) and iSchmider et all ([1991D sug- 
gest that the lo west-degree p- modes are the most likely 
to be detected. IGaulme et all (|201lD and iMarlevi (|1990f ) 
take the degree I = 8 to be an upper lim it and the Jo- 
vian observations of IGaulme et al.1 (|2011l ) are within the 
frequency regions [0.8, 2.1] mHz and [2.4, 3.4] mHz. In 
the specific case of Jupiter, and for I g [0, 8], these values 
loosely correspond to the ranges of radial order n G [0, 12] 
and [15,20], respectively. Theoretic ally, the asymptotic 
trends are m anifest for n > 4 or 5 ([Provost et al.1119931 ; 
IMarlevi ["19901. Thus, we focus on the p-modes and f- 
modes of low degree (I < 8) and relatively low radial 
order (n < 12). 

In §2, we summarize the theory of adiabatic nonradial 
oscillations of nonrotating spherically planetary models. 
We pres ent our nume r ical te chniq ue, closely based on the 
work of lUnno et al.l ([1989D and IChristensen-Dalsgaardl 
( 1997). To test the validity and precision of our code, we 
calculate the eigenfrequencies of f-modes, p-modes and 
g-modes of well-studie d polytrope models and compare 
our re sults to those of lChristensen-D alsgaar d~fe Mullanl 
(fl99l . 

In §3, we describe the giant planet models used in the 
article. We present our results in the case of Jupiter and 
Saturn, and compar e them to both obser yationnal data 
(IGaulme et al.ll201 If) and theoretical work (iProvost et al.l 
H993HMosser et al.lll994tlGudkova fc Zarkhovlll999ir ~We 
briefly focus on the dependence of the Jovian modal oscil- 
lations on the core mass of Jupiter and present the deriva- 
tives of the low-degree acoustic modes eigenfrequencies 
with respect to the core mass. 

In §4, we focus on the giant exoplanets. We present our 
results in terms of the characteristic frequency vq, the 
eigenfrequencies of low-degree f-modes, and the eigen- 
frequencies of low-degree p-modes across the giant exo- 
planet continuum. In separate subsections, we investi- 
gate their dependence on the planet mass, planet radius, 
and core mass. We also focus on the influence of a high 
fraction of heavy elements in the envelope, by using a 
high h elium mass fractio n, Y, as an approximate substi- 
tute ([Spiegel et al.ll2011l ). Finally, we briefly discuss the 
temporal evolution of the charatcristic frequency Uq for 
simple planetary models considered in isolation. 

2. METHODOLOGY AND TECHNIQUES 

2.1. Nonradial Oscillation Eigenvalue Problem 

We considered non-rotating, spherically-symmetrical 
planetary models. For adiabat ic, nonradia l oscil lations 
of such objects, it follows from lUnno et al.l ([19891 . chap. 
13) that the radial part of the displacement £ r , the Eu- 
lerian perturbations of the pressure p' , and the Eulerian 



perturbation of the gravitational potential $' take the 
form: 



f(t,r,6,<P) = f(r)Yr(e,< 



(1) 



for / = £ r , p' and A given oscillation mode is, thus, 
described by its azimuthal order m, degree radial order 
n. 

These variables are govern by a set of differential equa- 
tions and four boundar y conditions, two a t the surface, 
and two at the center ([Unno et al.l I1989T I . The corre- 
sponding set of equations is given in the first section of 
the Appendix. Since the azimuthal order m does not ap- 
pear in the governing equations, the eigenfrequencies are 
(2Z+l)-fold degenerate, and are fully described by their 
degree I and radial order n. 

This problem has to be numerically implemented to 
calculate the corresponding eigenfrequency v n _\ for a 
given mode. We detail the technique used in this paper 
in the Appendix, but summarize our overall methodology 
in the next subsections. 

2.2. Numerical Implementation 

Several n umerical techniqu es hav e been previously 
introduced dVorontsoy et all 119761 : lUnno et all 119891 : 
IChristensen-Dalsgaardl 119971 ). They are divided into 
shooting techniques and relaxation methods. In the 
shooting technique, solutions satisfying the boundary 
conditions are integrated separately from the inner and 
outer boundaries, and the eigenvalue is found by match- 
ing these solutions at an arbitrary interior fitting point. 
The second technique is to solve the equations together 
with the normalization condition, and all but one of the 
boundary conditions, using a relaxation technique; the 
eigenvalue is then found by requiring that the remain- 
ing boundary condition be satisfied. The shooting meth- 
ods are generally considerably faster than the relaxation 
techniques, but their precision decreases as th e degree 
I increases ([Christensen-Dalsgaardl 119971 120031 ). How- 
ever, since we consider only low-degree modes, a shooting 
method is quite suitable for our problem. Dimensionless 
variables are introduced (see the second section of the 
Appendix). In particular, the dimensionless frequency to 
is defined by: 



a 2 R? Att 2 v 2 R* 



GM 



GM 



(2) 



Solutions are obtained by integration using a fifth- 
order Runge-Kutta technique. To calculate the eigen- 
frequencies in a given freq uency range, we use a determi - 
nant method developed in[Christensen-Dalsgaard (|1997[ ) 
and fully described in the third section of the Appendix. 
Two linearly independent solutions are calculated from 
the center, and two from the surface. They are connected 
at an arbitrary inner boundary. The eigenvalues do not 
depend on its position. 

2.3. Mode Order 

As we calculate the eigenfrequencies and eigenfunc- 
tions of the acoustic modes for a given degree I, we 
determine their order n usi ng the following equation 
(|Christensen-Dalsgaardl[l997[ ): 



E 



n = > sign 

Xzl>0 



ill 



dyi\ 
dx I 



n 



(3) 
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where y\ and y^ are dimensionless variables denned by: 



2/1 = —, 
r 



V2- 



1 (t 



(4) 

where p is the density and g is the gravitational acceler- 
ation. 

In the definition of n, the sum is over the zeros x z \ in 
yi (excluding the center), where x = r/R is the relative 
radius. The value of no depends on the behavior of the 
solution close to the innermost boundary. If yi and j/2 
have the same sign at the innermost mesh point, exclud- 
ing the center, no = 0. Otherwise uq — 1. In particular, 
for a complete model that includes the center, as in our 
case, it follows from the boundary conditions at the cen- 
ter that riQ = 1 for radial oscillations and no = for 
non-radial oscillations. 

2.4. Results for a Polytrope Model 

In order to test the code, we compute the eigen- 
frequencies and eigenfunctions of polytropic mod- 
els. To compare our result s with the work of 
iChristensen-Dalsgaard fc Mullanl (|1994f ). we take the 
same radius and mass for our calculations: Rp = 
6.9599 x 10 10 cm and M P = 1.989 x 10 33 g. The eigenfre- 
quencies of f-modes, g-modes, and p-modes are given in 
T ables [H [H and El in p,Hz. Our frequencie s match those 
of lChristensen- Dalsgaa rd~fe Mullanl (|1994l ) to a precision 
of 10~ 5 or better, for all types of modes, for I £ [0, 3] and 
n £ [—20,25]. This gives us confidence in our calcula- 
tional method as we approach more complex models. 

3. JUPITER AND SATURN 
3.1. The Models 

The giant planet models we use for the calculations for 
Jupiter and Saturn consist of an adiabatic atmosphere, 
a hydrogen-helium envelope, and an olivine core. When 
a core is included, we explore < M core < 10 Mq for 
Jupiter and 9 < M core < 22 M® for Saturn. Both ranges 
are marginally consistent with the core accretion forma- 
tion models for these planets, which sugg est 10 — 20M ffi 
(ISaumon fc Guillotl l200l IPollackl [l996h . The hydro- 
gen/helium equ ation of state that we use for th is stud y 
is described in ISaumonT Chabrier. fc Van Horn! (fl995h . 
The transition between the atmosphere and the envelope 
has been smoothed to ensure the continuity of density, 
pressure, and sound speed. 

We build several models for both Jupiter and Saturn, 
using different core masses and helium fractions in the 
envelope. Table H] presents various parameters of these 
models: the helium mass fraction inside the envelope, Y, 
the mass of the core, M core (in Earth units), the cen- 
tral pressure p c (in Mbars), the central density p c (in 
cgs units) , and the characteris tic frequency vq (in pHz), 
defined in lProvost et al.l (|1993f ) using the following equa- 
tion: 



dr 

planet 

Cf) 



(5) 



static equilibrium. 

Figure [T] portrays the profiles of the density, p, the 
gravitational acceleration, <?, and the sound speed, cq, for 
models J4 and S2, defined in Table SJ Both density and 
sound speed are discontinuous at the core interface. In 
the case of Jupiter, for the models defined in TableHl the 
calculated characteristic frequencies, vq, are consistent 
with t he observational value measured bv lGaulme et al.1 
(|2TjTlT ): v Q = 155.3 ±2.2 p,Hz (see also Figure [6]). 

3.2. Oscillation Modes 

Figure [2] depicts the eigenfrequencies of models S2 and 
J4 for I £ [0, 3] and n up to 25. These results are given in 
the form of echelle diagrams based on the results of the 
asymptotic theo ry for low-deg r ee osc illations developed 
in, for example, iProvost et al.l (JT993). This theory pre- 
dicts that, for low degree I and large radial order n, the 
eigenfrequencies v n< i of p-modes are, to a first approx- 
imation, proportionnal to the characteristic frequency, 

v n ,i ^ (n + v . (6) 

An echelle diagram presents the ratio v n ,i/vo as a func- 
tion of the difference v n ^/v$ — (n + E[l/2]), where E is 
the floor function. Thus, it allows us to see the deviation 
from the approximate value. Both panels of Figure [2] are 
in qualitati ve agreement with previous numerical re sults 
for Jupiter (|Provost et al.lll993t IMosser et al.111994 and 
their Fig. 3d fc Fi g. 5a, respectively) and for Saturn 
(|Mosser et al.|[l994l their Fig. 5b). The Jovian periods 
of the acoustic fundamental tone and overtones with ra- 
dial order and degree up to 5 are given in Table [5] for 
the model J4. Though t he differences between our mod- 
els and those derived bv lGudkova fc Zarkhovl (|1999l ) pre- 
vent us from a precise comparison, it is clear that we find 
similar results to this previous work (see their Table 2). 

Figure [3] portrays the radial component of the eigendis- 
placement £ r for low-degree, lowest-order modal oscilla- 
tions of the model J4 of Jupiter, as a function of the rel- 
ative radius x = r/R. The radial displacement is shown 
for I = 0, 1, 2, and 5. It is taken equal to 1 m at the sur- 
face, for every mode. The behavior of the modes near the 
center is determined by the boundary condi tions of the 
speci fic numerical problem considered here ([Unno et al.l 
1989), which lead to the following relations, for r ~ 0: 



for I = 
for I > 1 . 



(7) 



where cq = U (dp a /dp Q ) a( i is the sound speed in hydro- 



Thus, for I = 1, the radial displacement does not neces- 
sary vanish near the center. In Figure [3l for the lowest 
degrees I, the presence of the dense core is directly visi- 
ble at the boundary with the envelope (x — 0.13, for this 
model). For higher degree (here, 1 = 5), the influence of 
the core is less obvious in the radial eigenfunctions, be- 
cause the amplitude of the radial displacement vanishes 
near the center, for every radial order n. 

The influence of the size of the core on the fre- 
quency spectrum of Jupiter has already been studied 
(|Provost et al.|[l993t iGudkova fc Zarkhovl [l999h . How- 
ever, no determination of the derivative of the eigen- 
frequencies with respect to the core mass has yet been 
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provided. Focusing on the f-modes and p-modes of 
Jupiter, we calculate their eigenfrequencies for various 
core masses, with the helium mass fraction fixed at 0.25 
in the envelope. An echelle diagram of the eigenfrequen- 
cies of Jupiter for I = 2 and for a few core masses is 
given in Figure 2] The spectra are very well separated 
for radial order n > 4, whic h indicates, as mentionned b y 
iVorontsov et all (|1989Q and lGudkova fc Zarkhovl (|1999t ). 
that the high-frequency acoustic oscillations of low de- 
gree I can be very useful in determining the structure 
and size of the core. 

We highlight the sensitivity to the core mass of the 
eigenfrequencies of low-degree acoustic oscillations. Fig- 
ure [5] shows the eigenfrequencies of such modes for 
Jupiter as a function of the mass of the core, M c , for 
n 6 [0, 10] and for I = 0, 1, 2 and 3. For every radial 
order n, the eigenfrequencies have been normalized by 
their coreless value: 

\ v n,l{M core ) 

VnJK-Mcore) — 7777, • (°) 

Vn,l{V>) 

This normalization allows us to compare the deviation 
of the frequencies from their coreless values as the core 
mass increases, regardless of their absolute value, which 
depends of the radial order n. 

For I = 0, as n increases the frequency becomes less 
sensitive to the size of the core. The normalized fre- 
quency /ii.o is by far the most affected by the variation 
of the core mass; its value decreases by more than 14% 
between the coreless version of Jupiter and the model 
with a 10- Earth mass core. For I > 1, the trend is very 
different; the derivative of the normalized frequency with 
respect to the core mass decreases from positive to neg- 
ative values as n increases. As a result, this derivative 
approximately vanishes for specific frequencies, for exam- 
ple /X3,i (FigEJ top right), ^2,2 (FigEJ bottom left), and 
A*2,3 (FiglU bottom right). All these frequencies vary by 
less than 0.5% over a core mass range from to 10 M®. 
For higher radial order, the influence of the core mass is 
more important; for example, ^10,2 varies by more than 
3.0% over the core mass range from to 10 M®. This 
increased sensitivity is consistent with the previous dis- 
cussion concerning the echelle diagram. The f-mode is 
also sensitive to the core mass, with a variation up to 
3.8% with core mass from to 10 M®, for I = 2 and 
3. We conclude that, for nonradial oscillations (I > 1), 
the sensitivity of the pulsation frequencies to the core 
mass is important both for the f-modes and for the high- 
radial order p-modes. For specific intermediate values 
of radial order n, the sensitivity vanishes over the whole 
core mass range shown. If detected, these modes would 
provide little evidence of the presence of a core. 

4. THE GIANT EXOPLANETS 

4.1. General description 

For a systematic look at the exoplanets cur- 
rently known, we u se the catalog developed by 
iSchneider et al.l (|2011[ ) and available at the URL 
http://www.exoplanet.eu. As of the 7 th of July, 2012, 777 
confirmed planets are listed in this catalog. We limit our 
study to the planets whose radius and mass have both 
been estimated. Furthermore, we focus on the giant exo- 
planet regime and, therefore, we select the planets whose 



radii are within the range 0.8 Rj < Rp < 2.1 Rj. In 
terms of mass, most of the detected giant exoplanets have 
masses less than 5 Mj, but the distribu tion has a long 
tail t owards masses larger than 10 Mj (jUdry fe Santosl 
|2007| ). Numerically, 88% of the selected exoplanets have 
masses less than or equal to 5 Mj, and 94% have masses 
less than or equal to 10 Mj. In the 10 to 20 Mj inter- 
val, it is difficult to fix a clear upper limit for giant ex- 
oplanets masses, because the plan et population and the 
brown dwarf popul ation overlap (jUdrv fc Santosl [20071 : 
iLeconte et all 12009( 1. Therefore, though we restrict the 
mass range studied, we are aware of the ambiguous sta- 
tus of the heaviest objects. This final set is composed of 
174 exoplanets. 

We calculate the characteristic frequency vq for each 
object of this group, using the techniques and models 
developed and described in sections 2 and 3. The he- 
lium mass fraction is fixed at 0.25 in the entire envelope 
and no core is added. The results are shown in Table 
[6l with the planets sorted from low to high mass. For 
the selected objects, the characteristic frequency range 
is 33 fj,Hz < vq < 815 [iHz. Around 1 Mj, vq is smaller 
than Jupiter's value for almost every object, since their 
radii are larger than 1 Rj. The spread of values is dra- 
matic around every mass within the range [0.17 Mj, 30 
Mj], which emphasizes the strong sensitivity of vq to 
planet parameters. 

To investigate the crossed dependence of Vq on the ra- 
dius and the planet mass, we calculate it for a wide range 
of radii and masses in the observed giant planet regime. 
Figure [5] portrays the corresponding results for 0.5 < 
M P < 10 Mj, and 0.95 < R P < 2.1 Rj. All the planet 
models are coreless, except for Rp = 1.0 Rj. For this 
radius value, we calculate the function i>o{Mp) for sev- 
eral core masses within the Jovian range < M core < 10 
M m . We place the ob served point for Jupiter, taken from 
iGaulme et"aH (|201lD . at M P = 1.0 Mj. As can be seen, 
our model is consistent with the observational data. For 
any fixed value of the planet radius, it is clear that Uq 
is an increasing function of the planet mass Mp. How- 
ever, the sensitivity of vq to the planet mass decreases 
as the planet radius increases. In order to quantify this, 
we fit the curves of Figure [5] with straight lines, in the 
high-mass regime (5 < Mp < 10 Mj), and calculate their 
derivatives. We find that, for a radius equal to 0.95 Rj, 
the corresponding derivative is 25 fj,Hz MJ 1 . At the ex- 
treme opposite end of the giant planet radius spectrum, 
for a radius equal to 2.1 Rj, the corresponding derivative 
is 8.0 fiHz MJ 1 . 

We investigate the separate influence of the radius, the 
mass, and the core mass on giant exoplanet pulsation 
modes. We focus on one parameter at a time. We se- 
lect three specific quantities to discuss this influence: the 
characteristic frequency the eigenfrequencies of low- 
degree f-modes, and the eigenfrequencies of low-degree 
p-modes. 

To determine the influence of planetary parameters on 
the characteristic frequency and on the frequency spec- 
trum of exoplanets, we calculate these for a wide range of 
each parameter (radius, mass, entropy, and core mass), 
all other things being equal. 
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4.2. Eigenfrequencies and characteristic frequencies of 
giant planets 

4.2.1. Dependence on the planet mass 

We build several planetary models with the radius fixed 
at Rp = 1.0 Rj and with various masses. We use coreless 
models, since we are here exploring the dependence on 
mass. Table [7] presents various parameters of these mod- 
els: the planet mass, Mp (in Jupiter units), the central 
pressure, p c (in Mbars), the central density, p c (in cgs 
units), the specific entropy, S (in kp /baryon), and the 
characteristic frequency, v$ (in fiHz). When the radius 
is fixed at 1.0 Rj, Uq is an increasing function of planet 
mass. We fit this function with a power law and obtain: 

(M \ ' 48 

v (M P ) ~ 164.0 x f piHz. (9) 

This dependence is due to the behavior of the the sound 
speed co(r), which is linked to the frequency vq by Equa- 
tion [5j Figure [7] depicts the profiles of the pressure po 
and the sound speed Co along the relative radius x = r/R, 
in hydrostatic equilibrium, for the first models of Table 
[71 As can be seen, at every level of the relative radius x, 
when the radius is fixed, both the pressure and the sound 
speed are increasing functions of the planet mass. Thus, 
given its definition (Eq. [5]), vq increases as the planet 
mass increases, all other things being equal. 

According to the asymptotic theory, we know that, for 
a given low value of the degree I, and for large radial 
orders n, vq is approximately the frequency gap between 
two modes of consecutive radial order: v n +i,l — v n,i ~ vq- 
Thus, at high frequency, the spacing between eigenfre- 
quencies increases when the planet mass increases, for a 
given value of the planet radius. Numerically, when the 
radius is fixed at Rp — 1.0 Rj, for an object of 1.0 Mj, 
we know that the high-frequency modes are separated 
by ^155 p,Hz, since this corresponds to the Jupiter case. 
For a 5.0-Mj planet, the frequency gap between high- 
frequency modes is ~350 p.Hz, and, at the end of the 
giant planet regime, for a planet mass of 15.0 Mj, the 
frequency gap exceeds 500 p,Hz. 

We now calculate the eigenfrequencies of the lowest- 
order p-modes for the objects defined in Table [7J and 
for I £ [0,8]. Figure [8] presents the corresponding eigen- 
values for n up to 12, as a function of the degree I, for 
objects with mass equal to 0.5, 1.0, 2.0, and 3.0 Mj, 
and a radius fixed at 1.0 Rj. The frequency spectra of 
the four planets appear more and more distinct from one 
another as we go up in frequency, and as we go up in de- 
gree /. Indeed, at low I, the f- modes of the four planets 
are close to one another, whereas the difference between 
the modes with the same n and I increases rapidly with 
frequency. Numerically, for I = 2, the f- modes of the 
four planets are all within the range 0.08 < ^0,2 _ 0.21 
mHz. For I = 2 and n = 12, the difference between the 
eigenvalues for Mp = 0.5 Mj and Mp = 3.0 Mj is more 
than 2.3 mHz. This discrepancy at high radial order n 
is, of course, due to the differences of the characteristic 
frequency vq, which is a measure of the frequency scale 
at low-degree I and high-order n. At high frequency, the 
value of £/q for each planet is clearly visible on Figure |U 

This increase of the frequency range continues as the 
planet mass increases beyond 3.0 Mj. To appreciate 



the difference numerically, we focus on I £ [0, 8] and 
n £ [0,7]. Figure [9] presents several low-order eigenval- 
ues, as a function of the planet mass, for various values 
of the degree I £ [0,8]. For the calculated modes, it 
appears from the calculations that the minimum in fre- 
quency is always obtained for (n, I) = (0, 2) (middle left 
panel), and the maximum is obtained for the highest n 
and I considered: (n,l) = (7,8) (bottom right panel). 
This statement is true for every value of the planet mass 
Mp in the selected range. On the bottom right panel 
(1 = 8), the functions 1*0,2 {Mp) have been added (black 
dashed line). Thus, the frequency range of the calculated 
modes is contained between the black dashed line (2*0,2) 
and the solid gold line, defined by (1*7,8 )■ Numerically, 
The low-degree, low-order eigenfrequencies of a 1.0- Mj 
planet are in the range [0.11,1.8] mHz, whereas the same 
eigenfrequencies of a 15- Mj object are in the range [0.50, 
6.4] mHz. 

4.2.2. Dependence on the planet radius 

We build several planetary models with the mass fixed 
at Mp = 1.0 Mj and with various radii. Table |8] presents 
various parameters of these models: the planet radius, 
Rp (in Jupiter units), the central pressure, p c (in Mbars), 
the central density, p c (in cgs units) , the specific entropy, 
S (in kp/baryon), and the characteristic frequency, i/q 
(in pHz). When the mass is fixed at 1.0 Mj, 2*0 is a 
decreasing function of the planet radius. We again fit 
this function with a power law and obtain: 

(R \~ 209 

2*0 (Ap) ~ 164.0 X f-^J 11Hz . (10) 

If we compare Equations 191 and [TU1 we see that, for the 
selected ranges of values, the dependence of vq on the ra- 
dius is significantly more important than the dependence 
on the mass. As explained in the previous subsection, vq 
is approximately the frequency gap between two modes 
of consecutive radial order, for a given low value of the 
degree I, and for large radial order n. Thus, this de- 
creasing behavior results in a diminution of the frequency 
gap between high-frequency modes. Numerically, we can 
see that this gap is around 40 fiHz for a 2.0-Rj planet 
(again, the mass is equal to 1.0 Mj), which is less than 
26% of the Jovian value. 

We calculate the eigenfrequencies of the lowest-order p- 
modes for the objects defined in TableEl and for I £ [0, 8]. 
Figure[l0]presents the corresponding eigenvalues for n up 
to 12, as a function of the degree I, for objects with a ra- 
dius equal to 1.0, 1.2, 1.4 and 1.6 Rj, and a mass equal to 
1.0 Mj. It appears that the remarks of the previous sub- 
section, which deals with the dependence on the planet 
mass, also apply to Figure [TUJ Indeed, the frequency 
spectra of the four planet models appear more and more 
distinct from one another as we go up in frequency, and 
as we go up in degree I . The frequency range and scale of 
the low-degree, low-order eigenvalues decrease with the 
planet radius, when the mass is fixed, whereas the same 
parameters increase with the planet mass, when the ra- 
dius is fixed. For instance, numerically, the low-degree, 
low-order eigenvalues of a 1.6-i?j planet are between 0.07 
mHz and 0.91 mHz, whereas, in the case of a 1.0-i?j 
planet, the same modes have eigenfrequencies between 
0.1 mHz and 2.6 mHz. 
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4.2.3. Dependence on the core mass 

Even in the cases of Jupiter and Saturn, the pres- 
ence and mass of a dens e core is still not proven. 
IGudkova fc Zarkhovl JT999) have shown that measure- 
ments of the pulsation modes of Jupiter could constrain 
the dimensions of the core. This is likely to be true for 
exoplanets, if and when their modes are measured. Many 
extrasolar gia nt planets appear smaller than the theor y 
would allow (jBurrows et al.l l2007t iGuillot et all 120061 ). 
This anomaly can be explained by the presence of heavy 
elements in a dense core, which shrinks the radii of these 
planets. One famous example is the case of HD149026b, 
whose measured radius and mass sugg est the presence of 
a core mass in the range 45 - 90 M e (|Sato et al.ll2005h . 

We calculate the characteristic frequency vq as a func- 
tion of the core mass for a selection of exoplanets for 
which the presence of a core has been inferred. The 
results are given in Figure II H which also includes the 
characteristic frequencies for Jupiter and Saturn, as 
a ref erence. The core mass ranges have been taken 
from iSaumon fc Guilloti (|2004l ) for Jupit er and Saturn 
ISato et all (120051) for HD149026b and IBurrows et all 
(|2007D for the other planets. It is clear that, in any 
case, with the radius and the mass of the object fixed, 
vq is a decreasing function of the core mass. This can 
be easily explained: the presence of a dense core reduces 
the sound speed in the center of the planet (see, for ex- 
ample, Figured]) which ultimately increases the integral 

§ planet c ano -: thus, diminishes vo, which is inversely pro- 
portionnal to the latter. However, the sensitivity to the 
mass of the core is not identical among the selected ob- 
jects. We can see that Saturn and HD149026b are much 
more influenced by the core mass than the others. This 
is due to their small radius and mass, compared to the 
other selected planets. Indeed, Saturn's radius is 0.83 
Rj, its mass is 0.30 Mj, HD149026b's radius is 0.72 Rj, 
its mass is 0.36 Mj whereas all the other planets have 
radii within the range [1.0,1.23] Rj and masses within 
the range [0.54,1.30] Mj. In this way, for planets with a 
small radius and mass, the determination of vq through 
observation can be a powerful tool to investigate the pres- 
ence of a dense core. For example, for our models of 
HD149026b, v$ loses more than 26% of its value between 
a 45-M® core model and a 90-M® core model. Thus, 
even a rough estimate of the value of vq might give us 
information on the core of this type of planet. 

To investigate the dependence of the low-degree, low- 
order eigenfrequencies on the core mass, we build several 
exoplanet models with the radius fixed at Rp = 1.0 Rj, 
the mass fixed at Mp = 1.0 Mj, and the core mass in the 
range 0-100 M ffi . Table IHpresents the various parameters 
of these models: the core mass, M c (in Earth units), the 
central pressure, p c (in Mbars), the central density, p c (in 
cgs units), the specific entropy, S (in kp /baryon), and 
the characteristic frequency, vq (in fiHz). The lowest- 
order eigenvalues of modal oscillations are given in Figure 
[T2l as a function of the degree I, for models with a core 
mass equal to 0, 10, 20, 30, 50, and 100 M®, and with 
the radius and mass fixed at the Jovian values. When 
the radius and mass are fixed, the eigenfrequencies are 
monotonic functions of the core mass, but their direction 
of variation depends on their degree I and their radial or- 
der n. For I > 2, the eigenvalues of the f-modes increase 



when the core mass increases. At high n, for every I, 
it can be seen that the eigenfrequencies are decreasing 
functions of the core mass, all things being equal. For in- 
stance, numerically, between the coreless model and the 
model with a core mass fixed at 100 M®, the eigenval- 
ues decrease by ~15% at high n (n S [10, 12]), for every 
/ 6 [0,8]. Thus, for a given degree I, the f-modes and 
the high-order p-modes, as functions of the core mass, 
have opposite directions of variation. Consequently, the 
frequency range for a given / shrinks when the core mass 
increases. If the low-degree, low-order modes were to be 
unambiguously identified by its spherical harmonic quan- 
tum numbers, the corresponding frequency range may 
constrain the presence of heavy elements in the deep in- 
terior. 

We have also constructed several exoplanet models 
with the mass fixed at Mp = 1.0 Mj, and the specific 
entropy fixed at S — 6.67 kp /bar yon, which is the spe- 
cific entropy of our coreless model with Rp = 1.0 Rj 
and Mp = 1.0 Mj. These models possess a core with a 
mass in the range 0-100 M s and this set of planet models 
approximately probes the situation in which, for a given 
planet mass and a given age, the presence of a dense core 
shrinks the radius. Table [TUl presents various parameters 
of these models: the core mass, M c (in Earth units), the 
planet radius, Rp (in Jupiter units), the central pres- 
sure, p c (in Mbars) , the central density, p c (in cgs units) , 
and the characteristic frequency, vq (in uHz). For the 
selected values of planet mass (Mp = 1.0 Mj) and en- 
tropy (S = 6.67 kp/baryon), the shrinking of the radius 
with core mass is visible. For instance, the planet radius 
decreases by ~20 % when a 100- M® core is added. The 
direct consequence of the reduction of the radius is the 
increase of the characteristic frequency v>q. Qualitatively, 
the increase of vq is consistent with its dependence on the 
planet radius, discussed in section 14.2.21 

For the fixed radius and entropy, the eigenvalues of 
the lowest-order modal oscillations are given in Figure 
[T2l as a function of the degree I, for models with a core 
mass equal to 0, 10, 20, 50, and 100 M®. When the 
mass and the entropy of the planet models are fixed, 
the frequency range of the low-degree, low-order p-modes 
decreases as the core mass increases. This is again a 
direct consequence of the shrinking of the planet radius 
as the core mass goes up, for given mass and entropy. 
However, as can be seen, the frequency spectra of the 
models defined by M c = M® and M c = 10 M® are 
quite similar for I S [0,3] and n S [0,12]. In particular, 
for / = 1 and 2, the eigenfrequencies of these two models 
differ by less than 1% for the frequency range considered. 
These similarities may be due to two opposite effects. 
First, we know that the presence of a 10- M® dense core 
shrinks the model radius by 2% for the selected values 
of planet mass and entropy (see Table [T0|) . When the 
planet mass is fixed, this decrease of the radius causes 
an increase of the low-degree, low-order eigenvalues, as 
already discussed in section 14.2.21 On the other hand, 
when the planet radius and the planet mass are both 
fixed, the presence of a dense core implies a decrease 
of the same eigenvalues. Figure [12] shows that, for the 
selected values, for low core mass and low degree I, the 
presence of a dense core compensates for the effect of the 
radius reduction on the eigenvalues. Nonetheless, it can 
be seen that, for higher degrees I (I > 4) and for higher 
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core masses (M c > 20 -M©), the effects of the radius 
reduction exceed the pure effect of the core. 

4.2.4. Dependence on the metallicity 

If not contained in a dense cor e, heavy elem ents can 
be laced throughout the envelope fGuillot 20Q5|) . In this 
spirit, we investigate the influence of heavy elements in 
the envelope itself, regardless of the presence of a core. 
Though there is no published robust equation of state 
that properly includes heavy elements beyond helium, 
we can mimic their presence by using a higher helium 
mass fraction than Y = 0.25 (|Guillotll2008t iSpiegel et all 
I20T1 . We assume the excess of helium mass fraction Ay 
(compared with the default value Yq = 0.25) is given by 
the value of the metallicity Z: 

AY ~ Z . (II) 

Using the value from lAsplund et al.l (|2009T ) . we take 
Zq = 0.014, to be the heavy element mass fraction in 
the Sun. In this way, an helium fraction of Y — 0.30 
mimics a metallicity equal to —3.6 x solar metallicity. 

We build three coreless models of exoplanets with a 
helium mass fraction of 0.25 and 0.30, respectively, with 
the radius and mass fixed at the Jovian values. The 
corresponding parameters, in particular the character- 
istic frequencies isq, are given in Table [TT] As can be 
seen, a higher helium mass fraction, hence a higher frac- 
tion of heavy elements in the envelope, tends to lower 
the value of Vo, all things being equal. Figure [TT] por- 
trays the low-degree, low-order eigenfreqdeuencies of the 
modal oscillations for the models of Table ITT1 It appears 
that the remarks made concerning Figure [T2l which deals 
with the dependence on the core mass, can be also made 
for Figure 1141 Such similarities are expected since, in 
both cases, it is the dependence on the global fraction of 
heavy elements that is under consideration. As can be 
seen on Figure [TTI for low I and very low n, the eigenfre- 
quencies slightly increase with the helium mass fraction 
y, whereas for higher I and n, the eigenvalues unam- 
biguously decrease with Y . As previously stated when 
discussing the presence of a dense core, if the low-degree, 
low-order modes were to be unambiguously identified, 
the corresponding frequency spectrum might give us a 
hint of the presence of heavy elements in the interior. 

4.2.5. Evolution of v§ with time 

To investigate the evolution of Uq for a given explanet, 
we build simple cvolutionnary models of exoplanets with 
a mass fixed at 0.5, 1.0, 2.0, 10, and 20 Mj. We use 
the default fo r malism and mod eling tools outlined in 
IBurrows et all ([19971 l200l [1991 . The planets are in iso- 
lation, which means that no stellar irradiation is taken 
into account. No core has been added. As the specific en- 
tropy decreases with time, the planet's radius decreases. 
Figure [15] portrays the evolution of vq and the planet 
radius up to 5 Gyrs for the five fixed planet masses con- 
sidered. The large early radii of the models result in 
small values for the characteristic frequencies, compared 
to their final values. Numerically, i/o(0) is between 16% 
and 22% of the final vq values, for the five models consid- 
ered here. As the radius stabilizes, so does v$. We fit the 
curves with straight lines in the region [2,5] Gyrs, and 
we calculate the corresponding derivatives. We find that 



the derivative increases with the planet mass, from —3.1 
LiHz Gyrs- 1 for M P = 0.5 Mj, to -8.2 fiHz Gyrs' 1 
for M P = 20 Mj. 



5. CONCLUSION 

We have calculated the eigenfrequencies and eigenfunc- 
tions of the pulsational modes of planets for a broad 
range of giant exoplanet models. In particular, we have 
investigated the dependence of the characteristic fre- 
quency vq, the eigenfrequencies of low degree f-modes, 
and the eigenfrequencies of low degree p-modes on sev- 
eral parameters: the planet mass, the planet radius, the 
core mass, and the helium fraction in the envelope. We 
provide the corresponding eigenvalues for a degree I up 
to 8 and a radial order n up to 12. We also present val- 
ues of v for 174 known giant exoplanets, and highlight 
the strong dependence on the radius, around any value 
of the planet mass. 

For Jupiter and Saturn, we find that our results are 
consi stent with both observationnal data (IGaulme et"aLl 
2TjTl and previous theoretical work (IProyost et al.lll993t 
Mosser et al.l [1991 IGudkova fc Zarkhovl 119991 7" In the 
specific case of Jupiter, we presented v$ and the low- 
degree, low-order eigenfrequencies of acoustic modes as 
a function of the core mass. We conclude that, for non- 
radial oscillations (I > 1), the sensitivity of the pulsation 
frequencies to the core mass is important both for the 
f-modes and for the high-order p-modes. However, for 
specific intermediate values of radial order n, this sen- 
sitivity is minimal over the whole range of core masses 
considered. 

Focusing on giant exoplanets, we find that the depen- 
dence of the characteristic frequency on the core mass is 
more important for small radii and masses. As an ex- 
ample, the characteristic frequency of HD149026b, with 
measured radius and mass of 0.72 Rj and 0.36 Mj, varies 
by more than 26% across the range of core masses con- 
sidered. We quantify the influence of the core mass on 
exoplanet models with arbitrary fixed mass and radius. 
For I E [0,8] and for frequencies up to 2.6 mHz, we find 
that eigenfrequencies shrink as the core mass increases, 
which is consistent with previous work on Jupiter. A big 
core (M c = 100 M®) induces a reduction in v n j of —15% 
for n > 10 compared to coreless values. We also develop 
an approach to quantify the influence on the eigenfre- 
quency spectrum of a high heavy-element fraction in the 
planet envelope. We find that, quantitatively, the pres- 
ence of heavy element in the envelope affects the eigen- 
value distribution in ways similar to the presence of a 
dense core. 

We have also quantified the influence of mass and ra- 
dius on the modal oscillations of giant exoplanets. For 
the selected values of I and n, we find that the pulsation 
eigenfrequencies depend strongly on both parameters, es- 
pecially at high frequency. This dependence can be mea- 
sured through 1/q. For the mass range 0.5 < Mp < 15 
Mj, and fixing the planet radius to its Jovian value, we 
find that i/ - 164.0 x (M P /Mj) 0A8 fiHz. For the ra- 
dius range from 0.9 to 2.0 Rj, and Axing the planet's 
mass to its Jovian value, we find that — 164.0 x 
(Rp/Rj) 2 09 nHz. These variations of vq directly af- 
fect the high-frequency spectrum of modal oscillations. 
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APPENDIX 

NONRADIAL OSCILLATION EIGENVALUE PROBLEM 

For adiabatic, non radial oscillations of non rotating spherically symmetrical planetary models, the governing differ- 
ential equations are ()Unno et al.lll989L eqs. 14.2 - 14.4): 

r^ dr c z \ a z ) pc z a^r^ 



p dr pc dr 



2 p' + {N'-^)Cr = -^r, (2) 



and 



1 d { 2 d&\ 1(1 + 1)^, _ ( p' N 2 



r 



r 2 dr \ dr J r 2 \ pc 2 ' g 

where c = (Tpo/ Po) 1 ^ 2 is the sound speed and the Lamb frequency, L\, is: 



$' = ^ Gp J- + — & , ( 3 ) 



and the Brunt- Vaisala frequency, N, is: 



Lf ^ (4) 



r = (dlnpo/dhi po)ad is the adiabatic exponent, G is the gravitation constant, g = GM r /r 2 is the gravitational 
acceleration, and r is the radius. Primed variables refer to the Eulerian perturbation at a given position; zero subscripts 
refer to the equilibrium value. 

Equations (Til) - (|3|) are the full fourth-order set of differential equations. There are four corresponding boundary 
conditions (|Unno etal.lll989L eqs. 14.8-14.11). At r = 0, 

Cr-^(^ + ^')=0 (6) 
a z r \ p J 



(7) 



and 

-j = ■ 

dr r 

At r — R, where R is the radius of the planet, 

— + = 8 

dr r 

and 

Sp = 0. (9) 

Here, S is the Lagrangian perturbation for a given fluid element. Equation ([9]) has various limiting forms, depending 
on the p hysical conditions at r — R. For the case in which the density and pressure vanish at the surface, ((9]) can be 
written (|Unno et al.lll989L eq. 14.12): 

£r - — = . (10) 

99 

This condition is valid whenever 

at r = R. Thus, we need to estimate the pressure scale height, H, for every input model of the planet. As an example, 
Saturn's value is about 40 km at 1 bar (Marley 1990), making this condition for a free boundary appropriate. The 
differential equations, boundary conditions, and a normalization condition at r — R, £, r /r = 1, comprise the eigenvalue 
problem. 
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EIGENVALUE CALCULATION: THE DIMENSIONLESS PROBLEM 

The eigenvalue differe ntial equations (E gsfT] - [3]) may be recast as four first-order differential equations with four 
dimensionless variables (jUnno et al.lll989l ). The variables are: 



2/1 : 

2/2 : 



r 

1 (t 

gr\p 

2/3 = —*', 
gr 

1 d& 
g dr ' 



2/4 = - 



Thc dimensionless variable 



r 

x = — 
R 



is used in place of r. The resulting four equations are as follows: 



2/2 + Vg?/3, 



= ( Cl ^ 2 - A*) yi + (A* - U + 1) y 2 + A*y 3 , 



a: -j— = (1 - t 7 ) 2/3 + 2/4, 
da; 



and 



^2/4 



= UA * Vi + + w + UV9] V3 + UyA - 



The dimensionless quantities are 



_ din M r 
s_ ~dTnT~ 



gr 



U 



1 dlnp Airpr 3 



Tdlnr M r 
d = (r/R) I (M r /M) 
, a 2 R 3 An 2 v 2 R i 



GM 



GM 



and 



A* = -rA = rg~ l N 2 . 



The dimensionless boundary conditions are 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 

(18) 

(19) 
(20) 
(21) 

(22) 



CiUJ 



2/i - 2/2 = at r = 



'2/3 ~ 2/4 = at r = 
(I + 1)2/3 + 7/4 = at r = -R 
2/1 - 2/2 + 2/3 = at r = R. 



(23) 

(24) 
(25) 
(26) 
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NUMERICAL IMPLEMENTATION 

Using the previous equations and boundary conditions, we build two linearly independent solutions that satisfy the 
appropriate boundary conditions at the center and two linearly independent solutions that satisfy the appropriate 
boundary conditions at the surface. These solutions are carried out using a fifth-order Runge-Kutta method with 
adjusted stepsize to ensure accuracy. The solution vectors y = (j/j) are, respectively: y c,1 (x), y c,2 (x), y s,1 (x), and 
y S ' 2 (a;), where the superscripts C and S denote a solution integrated from the center and the surface, respectively. 
A continuous match of the interior and exterior solutions at an arbitrary fitting point Xf requires the existence of 

non-zero constants iff' 1 , iff* 2 , K?' 1 , and iff' 1 , i G {1,2,3,4} such that (|Christensen-Dalsgaardll2003l) : 

K? A yf>\x f ) + K?> 2 y?< 2 (x f ) = K^yf'\x f ) + K?' 2 y?< 2 (x f ), (27) 

For all i 6 {1,2,3,4}. This set of equations has a solution only if the determinant 

A f (u 2 ) = \y c <\x f ) y c *(x f ) y s -\x f ) y s > 2 (x f )\ (28) 

vanishes. Hence, the eigenfrequencies are determined as the zeros of A/(w 2 ). This determinant has the advantage of 
behaving smoothly over the whole frequency spectrum and allows us to use a simple bisection method to find all the 
roots of Ay, while scanning a given interval of frequency. 

The roots of Af are supposed to be independent both of the choice of the fitting point and of the initial values of 
yC.i^ yC,2^ yS,i an( ^ yS> 2 . However, it is possible to control the amplitude of the determinant by using the regular 
solutions near the center and the surface, given in Unno et al. 1989 and characterized by: 

yi~x l ~ 2 forr-0, (29) 
yx~x~ l for r ~ R. (30) 
Thus, we define the following initial values for the center: 

C,X C,2 1-2 /oi\ 

Vi =Vx = x , (31) 

Vs ,1 = h -2/r\ (32) 

yf' 2 = /2 • (33) 

and for the surface: 

5.1 S.2 1-2 tnA\ 

Vi =Vx = x , (34) 
vl' 1 = 9i ■ vl' 1 , (35) 

5.2 S.l /on\ 

V 3 =92 ■ 2/i , (36) 

where /i, /2, gi and gi are arbitrary coefficients such that fx ^ f% and gx ^52- Then, for any particular frequency 
(actually w 2 ), the values of yf J and yf J for i = 2, 4 and j — 1,2 are fixed by the boundary conditions (|23|) - (|26|) . 
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TABLE 1 

ElGENFREQUENCIES OF P-MODES OF DEGREE AND 1 
FOR A POLYTROPE OF INDEX 3, IN llHz. THE RADIUS IS 

R P = 6.9599 x 10 10 cm and the mass is 
M P = 1.989 x 10 33 g. 



1 







1 




n 


1ms work 


CDM 


1ms work 


CDM 


1 


303.7754 


303.7755 


337.2152 


337.2152 


2 


411.5268 


411.5269 


463.5716 


463.5718 


3 


532.9179 


532.9181 


590.0692 


590.0694 


1 


658.1677 


658.1679 


716.6286 


716.6289 


5 


784.3664 


784.3667 


843.1062 


843.1066 


6 


910.7238 


910.7242 


969.4328 


969.4331 


7 


1036.9925 


1036.9929 


1095.5829 


1095.5832 


8 


1163.0917 


1163.0922 


1221.5526 


1221.5530 


9 


1289.0005 


1289.0010 


1347.3477 


1347.3484 


10 


1414.7209 


1414.7214 


1472.9793 


1472.9799 


11 


1540.2639 


1540.2645 


1598.4587 


1598.4594 


12 


1665.6437 


1665.6443 


1723.7983 


1723.7990 


13 


1790.8749 


1790.8756 


1849.0098 


1849.0105 


14 


1915.9720 


1915.9728 


1974.1042 


1974.1051 


15 


2040.9483 


2040.9489 


2099.0918 


2099.0928 


16 


2165.8149 


2165.8160 


2223.9823 


2223.9831 


17 


2290.5839 


2290.5848 


2348.7834 


2348.7844 


18 


2415.2641 


2415.2651 


2473.5035 


2473.5043 


19 


2539.8643 


2539.8653 


2598.1491 


2598.1500 


20 


2664.3924 


2664.3933 


2722.7266 


2722.7276 


21 


2788.8549 


2788.8558 


2847.2419 


2847.2430 


22 


2913.2577 


2913.2588 


2971.6997 


2971.7011 


23 


3037.6068 


3037.6078 


3096.1059 


3096.1067 


24 


3161.9068 


3161.9076 


3220.4627 


3220.4639 


25 


3286.1614 


3286.1626 


3344.7754 


3344.7767 


References. — 


- * Christciisen-Dalsgaard & Mullan 


41994 Table 2). 



TABLE 2 

ElGENFREQUENCIES OF P-MODES OF DEGREE 2 AND 3 
FOR A POLYTROPE OF INDEX 3, IN [iHz. THE RADIUS IS 
R P = 6.9599 X 10 10 cm AND THE MASS IS 

M P = 1.989 x 10 33 g. 



I 


2 




3 




n 


This work 


CDM* 


This work 


CDM* 


1 


390.1222 


390.1223 


428.8390 


428.8391 


2 


516.1991 


516.1992 


558.2978 


558.2981 


3 


643.0675 


643.0677 


686.8728 


686.8732 


1 


769.7799 


769.7802 


814.7025 


814.7027 


5 


896.2907 


896.2910 


942.0263 


942.0267 


6 


1022.6131 


1022.6135 


1068.9862 


1068.9866 


7 


1148.7618 


1148.7622 


1195.6634 


1195.6639 


8 


1274.7484 


1274.7491 


1322.1083 


1322.1090 


9 


1400.5845 


1400.5849 


1448.3544 


1448.3552 


10 


1526.2791 


1526.2795 


1574.4258 


1574.4265 


11 


1651.8419 


1651.8424 


1700.3401 


1700.3411 


12 


1777.2824 


1777.2825 


1826.1125 


1826.1136 


13 


1902.6072 


1902.6082 


1951.7557 


1951.7561 


14 


2027.8269 


2027.8277 


2077.2784 


2077.2793 


15 


2152.9475 


2152.9485 


2202.6912 


2202.6923 


16 


2277.9762 


2277.9776 


2328.0025 


2328.0032 


17 


2402.9209 


2402.9216 


2453.2187 


2453.2197 


18 


2527.7854 


2527.7867 


2578.3473 


2578.3482 


19 


2652.5775 


2652.5785 


2703.3935 


2703.3949 


20 


2777.3017 


2777.3023 


2828.3646 


2828.3655 


21 


2901.9614 


2901.9629 


2953.2632 


2953.2651 


22 


3026.5636 


3026.5648 


3078.0968 


3078.0983 


23 


3151.1104 


3151.1121 


3202.8689 


3202.8695 


24 


3275.6070 


3275.6085 


3327.5814 


3327.5828 


25 


3400.0566 


3400.0577 


3452.2396 


3452.2418 



Ref erences. — * IChristcnscn-Dalsgaar d & Mullanl 
(1994!, Table 2). 



Pulsation Frequencies and Modes of Giant Exoplanets 



TABLE 3 

ElGENFREQUENCIES OF G-MODES FOR A POLYTROPE 
OF INDEX 3, IN llHz. THE RADIUS IS 
R P = 6.9599 X 10 l0 cm AND THE MASS IS 
Mp = 1.989 x 10 33 g. 



I 2 3 



n 


This work 




This work 


ODM 


-20 


31.1798 


31.1797 


43.0341 


43.0356 


-19 


32.6437 


32.6439 


45.0003 


45.0009 


-18 


34.2535 


34.2535 


47.1556 


47.1552 


-17 


36.0311 


36.0313 


49.5267 


49.5272 


-16 


38.0054 


38.0053 


52.1528 


52.1519 


-15 


40.2096 


40.2099 


55.0737 


55.0718 


-14 


42.6886 


42.6884 


58.3407 


58.3400 


-13 


45.4951 


45.4953 


62.0230 


62.0230 


-12 


48.7014 


48.7010 


66.2057 


66.2053 


-11 


52.3970 


52.3972 


70.9964 


70.9960 


-10 


56.7061 


56.7065 


76.5390 


76.5389 


-9 


61.7959 


61.7959 


83.0266 


83.0267 


-8 


67.8987 


67.8991 


90.7245 


90.7243 


-7 


75.3534 


75.3535 


100.0041 


100.0058 


-6 


84.6648 


84.6649 


111.4182 


111.4181 


-5 


96.6263 


96.6264 


125.7916 


125.7915 


-4 


112.5543 


112.5543 


144.4538 


144.4530 


-3 


134.7963 


134.7963 


169.6560 


169.6563 


-2 


167.9278 


167.9279 


205.5287 


205.5286 


-1 


221.3677 


221.3677 


259.7578 


259.7578 



References. — * 

Christcnscn-Dalsgaard & Mullan ( 1994, Table 
4). 



TABLE 4 

Parameters of the models of Jupiter (J) and Saturn (S). 



Model 


Y 


Mcore (M©) 


p c (Mbars) 


p c (g cm 3 ) 


uq (pHz) 


Jl 


0.25 


0.0 


39.2 


4.14 


157.4 


.J 2 


0.30 


0.0 


41.7 


4.25 


153.0 


•J 3 


0.25 


5.0 


71.6 


19.4 


152.4 


J4 


0.25 


10.0 


96.1 


21.8 


151.2 


Si 


0.25 


13.0 


44.2 


16.1 


118.6 


S2 


0.25 


18.1 


59.8 


18.1 


115.0 


S3 


0.30 


18.1 


60.7 


18.2 


112.2 



TABLE 5 

Periods of p-modes for the J4 model (in min) . 



1 

n 





1 


2 


3 


4 


5 









138.38 


101.18 


84.49 


74.46 


1 


104.31 


59.44 


48.07 


41.52 


37.11 


33.93 


2 


45.53 


34.19 


30.13 


27.41 


25.42 


23.89 


3 


30.61 


24.93 


22.74 


21.21 


20.00 


19.00 


4 


23.50 


20.08 


18.55 


17.46 


16.60 


15.88 


5 


19.24 


16.92 


15.75 


14.94 


14.26 


13.69 
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TABLE 6 

Characteristic frequency v (in \xHz) for a set of 
coreless models using the estimated radius and mass 
(here in Jupiter units) of detected giant exoplanets. 



Planet 


ttplanet (ttj) 




uo (^Jrtz) 


Keplcr-9c 


0.823 


0.171 


94.40 


HAT- P- 18b 


0.995 


0.197 


65.01 


HAT-P-12b 


0.959 


0.211 


73.14 


Kcplcr-34b 


0.764 


0.220 


149.0 


WASP-29b 


0.792 


0.244 


141.9 


Kcplcr-9b 


0.842 


0.252 


114.1 


HAT-P-38b 


0.825 


0.267 


129.8 


WASP-39b 


1.270 


0.280 


47.87 


HAT- P- 19b 


1.132 


0.292 


61.16 


WASP-20b 


0.900 


0.300 


105.4 


WASP-21b 


1.210 


0.300 


54.29 


WASP-69b 


1.000 


0.300 


80.85 


HD-149026b 


0.718 


0.356 


229.4 


WASP-49b 


1.115 


0.378 


72.11 


WASP-63b 


1.430 


0.380 


44.52 


WASP-67b 


1.400 


0.420 


48.45 


Kcplcr-12b 


1.695 


0.431 


35.07 


Kcplcr-7b 


1.614 


0.433 


38.18 


WASP-llb 


1.045 


0.460 


92.36 


CoRoT-5b 


1.388 


0.467 


51.76 


WASP-31b 


1.537 


0.478 


43.59 


WASP- 13b 


1.365 


0.485 


54.39 


WASP- 17b 


1.991 


0.486 


33.51 


WASP-60b 


0.900 


0.500 


143.5 


WASP-42b 


1.080 


0.500 


89.56 


WASP-52b 


1.300 


0.500 


60.66 


WASP-6b 


1.224 


0.503 


68.68 


KOI-254b 


0.960 


0.505 


120.4 


HAT-P-lb 


1.217 


0.524 


70.94 


HAT-P-17b 


1.010 


0.534 


108.5 


CoRoT-16b 


1.170 


0.535 


77.76 


OGLE-TR-lllb 


1.077 


0.540 


93.89 


WASP- 15b 


1.428 


0.542 


52.87 


HAT-P-25b 


1.190 


0.567 


77.34 


WASP-62b 


1.390 


0.570 


56.96 


WASP-55b 


1.300 


0.570 


64.76 


WASP-25b 


1.260 


0.580 


69.56 


WASP-22b 


1.158 


0.588 


83.46 


WASP-34b 


1.220 


0.590 


74.97 


WASP-56b 


1.200 


0.600 


78.27 


WASP-54b 


1.400 


0.600 


57.66 


Kcplcr-8b 


1.419 


0.603 


56.37 


XO-2b 


0.973 


0.620 


130.3 


HAT-P-28b 


1.212 


0.626 


78.36 


Kepler- 15b 


0.960 


0.660 


139.9 


HAT-P-27b 


1.055 


0.660 


109.8 


Kcpler-6b 


1.323 


0.669 


67.80 


HAT-P-9b 


1.400 


0.670 


60.92 


OGLE-TR-lOb 


1.720 


0.680 


42.53 


HAT-P-4b 


1.270 


0.680 


74.19 


HAT-P-24b 


1.242 


0.685 


77.99 


WASP-59b 


0.900 


0.700 


181.6 


HAT-P-30b 


1.340 


0.711 


68.21 


HD-209458b 


1.380 


0.714 


64.62 


WASP-35b 


1.320 


0.720 


70.67 


CoRoT-4b 


1.190 


0.720 


87.62 


OGLE-TR-211b 


1.260 


0.750 


79.26 


TrES-1 


1.099 


0.761 


107.7 


HAT-P-33b 


1.827 


0.763 


40.7 


HAT-P-29b 


1.107 


0.778 


107.1 


WASP-68b 


0.900 


0.800 


196.7 


WASP-57b 


1.100 


0.800 


110.3 


CoRoT-9b 


1.050 


0.840 


126.6 


WASP-2b 


1.079 


0.847 


119.2 


HAT- P- 13b 


1.280 


0.850 


81.77 


WASP- 16b 


1.008 


0.855 


141.4 


WASP-lb 


1.484 


0.860 


61.93 


KOI-202b 


1.020 


0.880 


139.4 


WASP-23b 


0.962 


0.884 


163.6 


WASP-44b 


1.140 


0.889 


107.5 


WASP- 79b 


1.700 


0.890 


49.60 


XO-lb 


1.184 


0.900 


99.52 


WASP-28b 


1.120 


0.910 


113.3 


TrES-4 


1.706 


0.917 


50.07 



Note. — The helium mass fraction has been fixed at 0.25 in 
the entire envelope. No core has been added. Estimated radii 
and masses have been taken from http://www.exoplanet.eu 
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TABLE 6 

Characteristic frequency u (in \xHz) for a set of 
coreless models using the estimated radius and mass 
(here in Jupiter units) of detected giant exoplanets 



Planet 


K-planet \"-J 1 


M p l an et (™j) 


uq (^Jrtz) 


CoRoT-12b 


1.440 


0.917 


67.65 


WASP-41b 


1.210 


0.920 


95.95 


HAT-P-32b 


2.037 


0.941 


44.39 


WASP-7b 


1.330 


0.960 


80.61 


WASP-48b 


1.670 


0.980 


53.73 


WASP-45b 


1.160 


1.007 


110.4 


KOI-204b 


1.240 


1.020 


95.95 


WASP-26b 


1.281 


1.028 


90.00 


CoRoT-lb 


1.490 


1.030 


67.42 


WASP-24b 


1.104 


1.032 


125.3 


HAT-P-35b 


1.332 


1.054 


84.32 


HAT-P-6b 


1.330 


1.057 


84.70 


OGLE-TR-182b 


1.470 


1.060 


70.14 


HAT-P-5b 


1.252 


1.060 


95.89 


Qatar- lb 


1.164 


1.090 


114.2 


WASP-58b 


1.300 


1.100 


90.49 


CoRoT-19b 


1.450 


1.110 


73.68 


WASP-4b 


1.363 


1.121 


83.28 


HD-189733b 


1.138 


1.138 


122.9 


WASP-47b 


1.150 


1.140 


120.2 


WASP-78b 


1.750 


1.160 


54.10 


WASP- 19b 


1.386 


1.168 


82.39 


HAT- P- 3 7b 


1.178 


1.169 


115.3 


OGLE-TR-132b 


1.230 


1.170 


104.8 


OGLE-TR-113b 


1.110 


1.240 


136.2 


TrES-2 


1.169 


1.253 


121.6 


OGLE-TR-56b 


1.200 


1.300 


116.9 


CoRoT-13b 


0.885 


1.308 


562.3 


HAT-P-8b 


1.500 


1.340 


76.46 


WASP- 12b 


1.736 


1.404 


60.61 


WASP-50b 


1.153 


1.468 


136.3 


KELT-2Ab 


1.306 


1.486 


104.9 


WASP-65b 


1.300 


1.600 


110.1 


WASP-5b 


1.171 


1.637 


139.2 


WASP-37b 


1.136 


1.696 


151.9 


TrES-5 


1.209 


1.778 


135.3 


HAT-P-7b 


1.421 


1.800 


98.83 


HAT-P-36b 


1.264 


1.832 


125.2 


HATS- lb 


1.302 


1.855 


118.8 


TrES-3 


1.305 


1.910 


120.1 


HAT-P-15b 


1.072 


1.946 


187.1 


WASP-43b 


1.036 


2.034 


207.8 


WASP-61b 


1.240 


2.060 


138.6 


WASP-3b 


1.454 


2.060 


101.8 


HAT-P-23b 


1.368 


2.090 


115.0 


WASP-46b 


1.310 


2.101 


125.4 


Kcplcr-5b 


1.431 


2.114 


106.3 


HAT-P-22b 


1.080 


2.147 


193.1 


HAT-P-31b 


1.070 


2.171 


198.6 


HAT- P- 14b 


1.200 


2.200 


153.7 


KOI-428b 


1.170 


2.200 


162.5 


WASP-8b 


1.038 


2.244 


217.3 


CoRoT-21b 


1.300 


2.260 


132.4 


WASP-36b 


1.269 


2.279 


139.5 


WASP-66b 


1.390 


2.320 


118.0 


CoRoT-17b 


1.020 


2.450 


236.9 


Kepler- 17b 


1.312 


2.450 


135.8 


Qatar-2b 


1.144 


2.487 


182.0 


CoRoT-llb 


1.390 


2.490 


122.7 


WASP-53b 


1.200 


2.500 


164.3 


WASP-38b 


1.079 


2.712 


217.4 


CoRoT-lOb 


0.970 


2.750 


285.3 


CoRoT-23b 


1.050 


2.800 


235.6 


CoRoT-6b 


1.166 


2.960 


190.6 


WASP- 10b 


1.080 


3.060 


230.1 


HD-17156b 


1.095 


3.191 


227.5 


KOI-135b 


1.200 


3.230 


187.6 


CoRoT-2b 


1.465 


3.310 


129.9 


HAT-P-34b 


1.107 


3.328 


226.6 


CoRoT-18b 


1.310 


3.470 


163.9 


WASP-32b 


1.180 


3.600 


205.3 


HD-80606b 


0.921 


3.940 


380.3 


HAT-P-21b 


1.024 


4.063 


297.8 


HAT-P-16b 


1.289 


4.193 


186.6 



Note. — The helium mass fraction has been fixed at 0.25 in 
the entire envelope. No core has been added. Estimated radii 
and masses have been taken from http://www.exoplanet.eu 
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TABLE 6 

Characteristic frequency u (in \xHz) for. a set of 
coreless models using the estimated radius and mass 
(here in Jupiter units) of detected giant exoplanets 



ruanct 


ttplanet 


Mplanet \ M J ) 




OGLE2-TR-L9b 


1.614 


4.340 


126.5 


WASP-33b 


1.438 


4.590 


160.1 


HR-8799b 


1.100 


7.000 


328.2 


WASP- 14b 


1.281 


7.341 


252.4 


CoRoT-14b 


1.090 


7.600 


347.5 


KOI- 13b 


1.830 


8.300 


158.2 


Kepler-14b 


1.136 


8.400 


336.8 


HAT-P-2b 


0.951 


8.740 


490.3 


Keplcr-30c 


1.290 


9.100 


278.3 


SWEEPS- 11 


1.130 


9.700 


364.8 


HR-8799c 


1.300 


10.000 


288.1 


HR-8799d 


1.200 


10.000 


332.0 


WASP- 18b 


1.165 


10.430 


357.5 


lRXS1609b 


1.700 


14.000 


216.3 


HN-Pegb 


1.100 


16.000 


487.6 


Kcplcr-30d 


0.960 


17.000 


639.1 


KOI-423b 


1.220 


18.000 


432.3 


2M-2140+16b 


0.920 


20.000 


742.3 


GQ-Lupb 


1.800 


21.500 


264.3 


CoRoT-3b 


1.010 


21.660 


653.3 


2M-2206-20b 


1.300 


30.000 


501.3 


2M-0746+20b 


0.970 


30.000 


815.6 



Note. — The helium mass fraction has been fixed at 0.25 in 
the entire envelope. No core has been added. Estimated radii 
and masses have been taken from http://www.exoplanet.eu 
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TABLE 7 

Parameters of coreless exoplanet models of various masses with a radius 
fixed at Rp = 1.0 Rj. The helium mass fraction in the envelope is 0.25. 



M P (Mj) 


p c (Mbar) 


Pc 


(g cm 


3 ) Specific entropy S (kp /baryon) uq (p-Hz) 


0.5 


9.8 




2.0 


6.9 107.3 


1.0 


39.2 




4.1 


6.7 157.4 


2.0 


162.5 




8.6 


6.5 225.8 


3.0 


383.8 




13.6 


6.6 274.3 


5.0 


1176 




24.7 


7.1 345.4 


10.0 


5612 




55.8 


8.2 467.2 


15.0 


14075 




89.1 


8.9 560.5 










TABLE 8 


Parameters of coreless exoplanet models for various radii with a mass fixed 


ar Mp = 1.0 Mj. 


The helium mass fraction in the envelope is 0.25. 


Rp (Rj) 


p c (Mbar) 


Pc 


(g cm 


~ [i ) Specific entropy S (kp /baryon) vq (fJ.Hz) 


0.9 


48.2 




5.0 


3.5 220.2 


1.0 


39.2 




4.1 


6.7 157.4 


1.2 


28.7 




3.0 


8.7 102.0 


1.4 


22.5 




2.3 


9.7 74.57 


1.6 


17.9 




1.8 


10.3 58.47 


1.8 


14.2 




1.5 


10.7 47.78 


2.0 


11.5 




1.2 


10.9 40.08 










TABLE 9 


Parameters of exoplanet models for various core masses M c with a planet 


MASS FIXED AT Mp = 


1.0 Mj and a planet radius fixed at Rp = 1.0 Rj. The 




HELIUM MASS FRACTION IN THE ENVELOPE IS 0.25. 


M c (M®) 


p c (Mbar) 


Pc 


(g cm 


3 ) Specific entropy S (kp /baryon) uq (fJ.Hz) 





39.2 




4.1 


6.7 157.4 


10 


96.9 




21.8 


7.0 151.2 


20 


145.7 




25.8 


7.3 149.4 


30 


198.9 




29.2 


7.6 147.5 


50 


307.9 




35.5 


8.1 143.0 


100 


645.4 




49.5 


9.2 134.1 










TABLE 10 




Parameters of exoplanet models for various core masses 




M c WITH 


A planet mass fixed AT Mp = 1.0 Mj and a 



specific entropy fixed to 6.67 kp /baryon (value for a 

CORELESS MODEL WITH Rp = 1.0 Rj AND Mp = 1.0 Mj). THE 
HELIUM MASS FRACTION IN THE ENVELOPE IS 0.25. 



M c (M e ) 


Rp (Rj) 


p c (Mbar) 


p c (g cm 3 ) 


u (pHz) 





1.0 


39.2 


4.1 


157.4 


10 


0.98 


100.1 


23.0 


160.0 


20 


0.96 


152.2 


27.3 


166.8 


30 


0.94 


209.1 


31.2 


173.0 


50 


0.90 


329.5 


38.1 


187.3 


100 


0.81 


680.6 


53.0 


223.6 



TABLE 11 

Parameters of coreless exoplanet models for some helium mass 
fraction y in the envelope with a planet mass fixed at mp = 1.0 mj 
and a planet radius fixed at rp = 1.0 rj. 



Y 


Pc (Mbar) 


p c (g cm 3 ) Specific entropy S (kp /baryon) 


u (liHz) 


0.25 
0.30 


39.2 
41.7 


4.1 6.7 
4.3 6.9 


157.4 
153.12 
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Fig. 1. — Left panel: Distribution of density po (g cm 3 ), gravitational acceleration g (cm s 2 ), pressure po (Mbar), and sound speed 
co (km s —1 ) as a function of the relative radius r/R for the model J4 of Jupiter. Right panel: The same in the case of the model S2 of 
Saturn. One has to notice that the scales of the Y-axis are different on the two figures. 




.0.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 



v„ , / v - (« - E[l/2]) v„ , / V - (n - E[l/2]) 

Fig. 2. — Left panel: Echelle diagrams of the eigenfrequencies of Jupiter (model J4) for £ £ [0,3] and n £ [0,25]. The characteristic 
frequency, vq, is 151.2 /iHz (1 = 0: +, I = 1: X, I = 2: *, I = 3: □). Right panel: The same in the case of the model S2 of Saturn. The 
characteristic frequency, uq, is 115.0 /iHz. 
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Fig. 3. — Radial component of the eigendisplacement £ r for low-degree, lowest-order modal oscillations of the model J4 of Jupiter, which 
has a core mass equal to M c = lOiUg . The radial displacement is taken equal to 1 m at the surface; in other words, the radial displacement 
£ r is normalized by its value £i? at the surface. Top left: Z = 0; top right: 1=1; bottom left: 1 = 2; bottom right: / = 5. 
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-0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

v„,,/v -(«-£[Z/2]) 

Fig. 4. — Echelle diagrams of the Jovian eigenfrcqucncics calculated for different core masses with uo = 155 fiHz and 1 = 2 (In Earth 
units, the corresponding core masses are M core = 0: +, M core = 5: X, M core = 7: *, M core = 10: □). 
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Jupiter, / = 1 




Jupiter, 1 = 2 




Jupiter, 1 = 3 



AM®) 




Fig. 5. — Top left panel: Eigcnfrcquencies of p-modes for Jupiter as a function of the core mass for / = 0. For every radial order n, the 
cigenfrequencies have been normalized by their coreless value: fi n i(M core ) = v n i(M core )/i/ n ;(0). We assume that < M core < lOAf^. 
Top right panel: same as for the top left panel, but for 1 = 1. Bottom left panel: same as for the top left panel, but for / = 2. Bottom right 
panel: same as for the top left panel, but for I = 3. 
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I i i i i i i i i i I 

1 23456789 10 

M Pkmet (Mj) 



Fig. 6. — The charateristic frequency, uo, as a function of the planet mass, for various planet radii. The helium mass fraction in the 
envelope has been set to 0.25. The models are all coreless except for Rp = 1.0 Rj. For this radius value, we calculate the function vo(Mp) 
for various core masses between (solid red line) and 10 M(g (solid purple line). The figure includes an insert which zooms into the 
range 0.9 < Mp < 1.1 Mj. The gray area in the insert depicts the frequency range reached by uo for the models with Rp = 1.0 Rj and 
< M core < 10 Mg). The observed point for Jupiter has been added, both on the figure ( black cross) and on the insert (with errorbars). 
The value of its characteristic frequency is taken from the measurements oflGaulmc ct al. ( 201l|). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



r/R r/R 

Fig. 7. — Left panel: distribution of sound speed co (km s _1 ) as a function of the relative radius r/R for coreless exoplanet models 
with various masses, with the radius fixed at Rp = 1.0 Rj. Right panel: distribution of pressure po (Mbar) for the same models, with 
a logarithmic scale for the Y-axis. At every radius r, both the sound speed and the pressure are increasing functions of the planet mass, 
when the planet radius is fixed at Rp = 1.0 Rj. 
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Fig. 8. — Low-order eigenfrequencies of oscillation modes (in mHz) as a function of the degree I for coreless exoplanet models with various 
masses, with the radius fixed at Rp = 1.0 Rj. The helium mass fraction in the envelope is 0.25. For the sake of clarity, the eigenvalue 
dots of the different models has been seperated and drawn on both sides of each integer degree I. 
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Fig. 9. — Eigenfrcqucncics of oscillation modes (in mHz) as a function of the planet mass for coreless exoplanet models with the radius 
fixed at Rp = 1.0 Rj, and for various values of the degree I. For the represented modes, for each value of the planet mass Mp, the 
frequency minimum and maximum are obtained for (n,l) = (0,2) (middle left panel) and (n,l) = (7,8) (bottom right panel), respectively. 
On the bottom right panel (I = 8), the functions vo 2\Mp) has been added (black dashed line). Thus, the frequency range of the calculated 
modes is contained within the solid gold line, defined by (n, /) = (7, 8), and the black dashed line, defined by (n, /) = (0, 2). 
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Fig. 11. — The charateristic frequency vq as a function of the core mass, M c , for various exoplanets. The helium mass f raction has been set 
equal to 0.25 in the e nvelope. The estim ated range of core masses have been taken from various sources; Jupiter, Saturn: Saumon & Guillot 
pjQlT) ; HD149026b: ISato et al.l J2005T) ; for all the other planets: IBurrows et~aTI (f2007Tl . 
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Fig. 12. — Low-order eigenfrequencies of oscillation modes (in mHz), as a function of the degree I, for exoplanet models with various 
core masses, with the radius fixed at Rp = 1.0 Rj and the mass fixed at Mp = 1.0 Mj. The helium mass fraction in the envelope is 0.25. 
For the sake of clarity, the eigenvalue dots of the different models has been seperated and drawn on both sides of each integer degree I. 
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Fig. 13. — Low-order eigenfrequencies of oscillation modes (in mHz), as a function of the degree I, for exoplanet models with various 
core masses, with the radius fixed at Rp = 1.0 Rj and the specific entropy fixed at S = 6.67kg /baryon. The helium mass fraction in the 
envelope is 0.25. For the sake of clarity, the eigenvalue dots of the different models has been seperated and drawn on both sides of each 
integer degree I. 
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Fig. 14. — Low-order eigenfrequencies of oscillation modes (in mHz), as a function of the degree i, for exoplanet models with two different 
helium mass fractions Y in the envelope, with the radius fixed at Rp = 1.0 Rj and the mass fixed at Mp = 1.0 Mj. For the sake of clarity, 
the eigenvalue dots of the different models has been seperated and drawn on both sides of each integer degree I. 
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Fig. 15. — Left panel : the charateristic frequency vq as a function of time for various coreless planetary models, characterized by their 
constant mass (Mp = 0.5, 1.0, 2.0, 10.0 and 20.0Afj). The helium mass fraction has been set to 0.25 in the envelope. The planets are 
considered in isolation during their evolution, which means that no irradiation is taken into account. Right panel: The corresponding 
evolution of the planet radius for the same coreless planetary models, in Jupiter units. The X-axis does not begin at the origin. 



